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Abstract 

This paper establishes the universality of parametric correlations of eigen- 
functions in chaotic and weakly disordered systems. We demonstrate this 
universality in the framework of the gaussian random matrix process and 
obtain predictions for a number of parametric correlators, one of them ana- 
lytically. We present numerical evidence from different models that verifies 
our predictions. 
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The statistical fluctuations of spectra and wave functions in complex systems are well 
known to conform to the predictions of random matrix theory (RMT) These systems, 
whose common feature is their non-integrability, range from single-particle systems exhibit- 
ing chaotic behavior |2[ , such as ballistic quantum dots with an irregular confining potential 
||, to interacting many-particle systems, such as strongly correlated electron models || and 
atomic nuclei [f|. RMT predictions hold also for electron systems with a random impurity 
potential which is sufficiently weak to allow for diffusion. According to RMT, the eigenfunc- 
tions at a given point in space are gaussian random variables whereas if the energy 

levels are scaled by the mean spacing, the spectral correlations become independent of the 
details of the system and obey the Wigner-Dyson spacing distribution. 

It has recently been discovered that when these systems are allowed to depend on a 
parameter (e.g. an external field), the correlations between spectra belonging to different 
parameter values become universal upon an appropriate scaling of the parameter |7|], ||, 
||. The scaling factor turned out to be the RMS of the level velocity divided by the mean 
spacing. 

The purpose of this paper is to establish the universality of parametric correlations of the 
eigenfunctions in these systems. This universality is deduced within a general framework 
which we find suitable for this discussion: the gaussian process (GP) [|K|, a random matrix 
process corresponding to each one of Dyson's three gaussian ensembles (GE). We then 
concentrate on the case of conserved time-reversal symmetry and present two parametric 
correlators involving the eigenfunctions. We demonstrate the universality of these correlators 
by comparing their RMT predictions, obtained from an appropriate GP, with results of 
numerical simulations for a chaotic system and a disordered system. Finally we introduce 
a third correlator, which we are able to calculate analytically, and thus provide an explicit 
example to the general discussion of scaling. The universality in this case is again verified 
by numerical simulations. 

Dyson [JTTJ] showed that there are only three possible types of gaussian ensembles that can 
describe a physical system, depending on its symmetry: orthogonal (GOE), unitary (GUE) 
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and symplectic (GSE). When the system depends on a parameter x, and respects the same 



symmetry for all values of x, it has been proposed [JTOj] that its statistical properties may 
be described by corresponding gaussian processes termed respectively gaussian orthogonal 
process (GOP), gaussian unitary process (GUP), and gaussian symplectic process (GSP). 
A GP is a set of random N x N matrices H(x) whose elements are distributed at each x 
according to the appropriate GE with a prescribed correlation among elements at different 
values of x: 



Hij(x) = 



H tJ {x)H kl {x') = ^f{x,x')g%, (1) 

where = dikdji + Sufyk an d 9ifkP = ^u^jk- We will be concerned with stationary 

processes for which the process correlation function f(x,x') = f(\ x — x' |) and is normalized 
such that /(0) = 1. A GP is completely determined by its first two moments defined in 
([!]) and has the useful property that the joint probability distribution of any finite number 
of matrices H(x), H(x'), H(x"), ... is gaussian. In particular, at any x we have P [H(x)] oc 
exp [— /3TrH(x) 2 /2a 2 ]. The case (3 = 1 corresponds to conserved time-reversal symmetry 
and H{x) are real symmetric, whereas for (3 = 2 this symmetry is broken and H(x) are 
complex hermitean. As usual in applications of RMT we are interested in the limit N — > oo. 

An example of a gaussian process is Dyson's Brownian motion model [12|] where f(x) = 
exp(— 7 | a; |). This model has been used to relate the above-mentioned spectral correlators 
to the spacetime correlations of the Sutherland-Calogero-Moser system 0,|ll]. 



Considering now an arbitrary two-point correlation function c(x, x') = 0(x)0(x') where 
0(x) is some observable involving the spectrum and eigenfunctions of H(x), it can be shown 
p~0j , [15[1 that c(x — x')/c(0) depends, apart from on a, on the combination N[l — f(x — x')]. To 
make a correspondence with a particular system we note that a sets the mean level spacing 
A near the center of the spectrum through a/A = y2N /n. Absorbing a into H(x) in ([I]) 
is thus achieved by scaling the energies Ei by A to get the "unfolded" energies = Ei/A, 
leaving c(x — x')/c(0) independent of a. Next, for x' near x we expand 



f(x - x) ss 1 - k(x - x'f (2) 



and, following ||, we consider the variance of the level velocity (dei/dx) 2 . Writing 
dEi/dx ~ (ipi(x) | H(x') — H(x) \ ipi(x)) j[x' — x), where | ipi(x)) are the eigenstates of 
H(x), we calculate its variance using the joint two-matrix distribution P [H(x), H(x')]. The 
calculation is performed in two steps. First, using the conditional distribution for H(x') 
given H(x) which can be shown to be 

P [H(x') | H(x)} = P [H(x), H(x')\ jP [H(x)\ 

ex exp {-/3Tr [H(x') - fH(x)] 2 /2a 2 (l - f)} (3) 
with / = f(x — x'), we average over H(x') keeping H(x) fixed to get 



'dE,\ 2 1 



dx I (x — x 



l\2 



^(l-/ 2 ) + (l-/) 2 ^) 2 



(4) 



Second, taking the limit x' —>■ x and using (0), the second term on the r.h.s. of (^) vanishes 
and we obtain the following expression for the non-universal quantity k: 



Thus after the scaling 



x — ► x 



1/2 



(dei/dx) 2 x (6) 



we get 



and all correlators, being determined by the combination N(l — f), become universal as 
functions of x — x'. 

A few remarks are in order. First, the scaling (|6|) is identical to that found in || for 
spectral correlators. Second, this scaling was derived here under the assumption that the 
second derivative of / is the first non- vanishing one (see (0)). This is not always the case 



as exemplified by Dyson's Brownian motion model for which / « 1 — 7 | x |. The more 
general case is discussed elsewhere ||10| . Third, the form ([!]) of / implies that the typical 
correlation length x scales like ^= in the GP, as was shown in jl4j] and will be seen explicitly 
in the analytical expression for the correlator 0(0;, x — x') derived below. Fourth, and most 
importantly, we note that previous treatments of parametric correlations had to invoke the 



supersymmetry method || or Dyson's Brownian motion model |13| and could demonstrate 
universality of spectral correlators only (however, see 0]). In the GP framework, on the 
other hand, the universality of all correlators emerges quite simply. In particular, correlators 
involving the eigenfunctions are universal. 

We shall demonstrate the last point by investigating two such quantities, the averaged 
parametric overlap 



o(x - x ) = \ (i[>i(x) I ipi(x')) | 2 (8) 
and the projection correlator 



p(x-x') = (<p\^(x))(<P\Mx')) ■ (9) 

Both measure the decorrelation of wave functions as their separation along x increases. 
o(x — x') gives the overlap of wavefunctions at different x, whereas p(x — x') provides the 
correlation between their components along a fixed normalized vector | <ft). To see that p(x — 
x') is independent of the choice of | 0), notice that had we chosen instead | <fi') — U \ (j>) for 
some unitary U, we could have rotated H(x) by W at each x without affecting the probability 
measure (see ([|)), thereby recovering the original result for p(x — x'). In particular, for 
0) =| r) the projection correlator describes the correlation between eigenfunctions belonging 



to different values of a; at a given space point r: p(x — x') = ipf (r)ifjf (r) . For the orthogonal 
case ipi(r) are real and determined up to a sign, which is easy to keep fixed as x varies. 

We cannot provide analytical expressions for these correlators. However, since they are 
universal we can use any gaussian process to obtain a theoretical prediction for them from 
a random matrix simulation. A simple GOP is generated by H(x) = ificosx + i^sinx, 



where H\,H2 are independent GOE matrices. This process is stationary with f(x — x') = 
cos (x — x') and the scaling (§) is readily found to be x — > x = (y2N /tt)x. The theoretical 
curves, generated by a simulation of random matrices with N = 150, are given by the dashed 
lines in Fig. |1| 

To verify the universality we studied o(x — x') and p(x — x') in both a chaotic system and 
a disordered system. The first is the interacting boson model (IBM), a many-body system 
used to describe collective states of medium and heavy mass nuclei [l7fl . Its constituents are 



bosons which model nucleon pairs coupled to angular momentum of or 2. Depending on 
two parameters Xi V t ne IBM can be integrable or non-integrable. The time-dependent mean 
field equations, obtained in the limit of an infinite number of bosons and which constitute 



the classical limit, are correspondingly regular or chaotic JT^]. We have calculated the above 
correlators as a function of x m the regime 77 = 0, —0.8 < x < —0.5 where the system is 
almost fully chaotic. Since the total spin J of the nucleus is a conserved quantum number, 
we can study the correlations for different values of the spin. Results for J = 2 and J = 6 
with 25 bosons are displayed in the right panel of Fig. [I]. 

The second system we studied was the Anderson model, a two-dimensional lattice Hamil- 
tonian with on-site disorder and nearest-neighbor hopping. The site energies W%, measured 
in units of the hopping term, are uniformly distributed in [— W/2,W/2] where W controls 
the transition between the diffusive and localized regimes. We considered the cases of cylin- 
drical geometry with W = 4 and toroidal geometry with W = 2, introducing a parametric 
dependence by adding a step potential of strength x along one of the lattice directions [§ . 
Results for a 27 x 27 lattice are presented on the left panel of Fig. |I[ Both correlators in 
both systems are in excellent agreement with the GOP prediction. 

The last part of this paper introduces a third correlator which we are able to calculate 



analytically using Efetov's supersymmetry method JT9[. This correlator, related to the 



parametric overlap (M) but involving both energies and eigenf unctions, is given by 



o{uj,x — x) = — = . (10) 

Eij<5(e*W + u-e j (x')) 



It measures the averaged parametric overlap of eigenfunctions whose corresponding energies 
are separated by u in units of the mean spacing. We now outline the calculation of the 
numerator of (|10|), denoted o(u>, x — x') (details will be published elsewhere |15[). In order to 



employ the supersymmetry method we first express it in terms of Green functions: o(u, x — 

where e = e ± i5. Each 



x') = Tr^ReTr 



G(e-,x)G(e+ + u, x') - G{e-,x)G{e~ + u, x') 
Green function is then written as an integral over 4A r -dimensional graded vector \I> '. After 
performing the GOP averaging, a Hubbard- Stratonovich transformation and a subsequent 
integration over the ^-variables result in an integral over a 16-dimensional graded matrix 
R. This integration is carried out in the saddle-point approximation, which is exact in the 
infinite- N limit. In this limit the quadratic corrections decouple and we are left with an 
integral over the saddle-point manifold: 

o{u, x-x') = ^Re Jd[Q]P [Q] exp (F [Q]) (11) 

where Q is a 8 x 8 graded matrix and 

P[Q] = Trg(QSl*)Trg(QS 2 b Z) , 

F[Q] =t^(u + 2t5)Tig(QA) + ^NK(x-x') 2 TTg[Q,A] 2 . (12) 



Using Efetov's parametrization of Q [[0| most of the integrals are elementary and the result 
is given by 

o{u, x-x') = ^Re JD[X]P [A] exp(F [A]) , (13) 

where 

r« roo rl (1 - \ 2 )d\ 1 d\ 2 d\ 



/POO POO fl 
/ /_ 



h h J-i (A? + A| + A 2 -2A!A 2 A-1) 2 ' 
P [A] = 2XlX 2 2 - Xl - X 2 2 - A 2 + 1 , 

F [A] =tjr(uj + 2%8) (AiA 2 - A) - Nk{x - x') 2 (2A 2 A 2 - X\ - A 2 - A 2 + 1) . (14) 



Note the combination N(x — x') 2 in o(cu, x — x') which explicitly suggests the yN scaling of x 
mentioned above. The scaling @ removes the dependence on the system-specific k through 



the substitution Nk(x ) 2 . Our result for the denominator of (|10|) , after 

scaling, becomes identical to the level density correlator calculated in [f|. An expression 
for o(u, x — x') can similarly be derived for the case of broken time-reversal symmetry by 
GUP-averaging the Green function product ||15|| . 

We verified our derivation of o(uj, x — x') by a comparison with a GOP simulation and 
confirmed the universality of this correlator by studying it in the above two cases of the 
Anderson model. We remark that the expression in (Ti3|)-(|T4|) corresponds to a regularization 
of (|1(]) by convoluting both the numerator and the denominator with a Lorentzian of width 
5. The same operation, which amounts to smearing the 5-functions, was performed in the 
numerical computations. The results are displayed in Fig. for 5 = 0.05 in units of the mean 
spacing and various values of uj and are in excellent agreement with the GOP prediction. 

In conclusion, we have shown that the concept of the gaussian process is particularly suit- 
able for a discussion of parametric correlations in chaotic and disordered systems and for 
demonstrating their universality. We established the existence of universal parametric cor- 
relations of eigenfunctions in these systems, provided predictions for three such correlators, 
and verified them in different models. 

This work was supported in part by the Department of Energy Grant DE-FG02- 
91ER40608. 
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FIGURES 

FIG. 1. The eigenfunction correlators o(x — x') (top) and p(x — x') (bottom) as a function of 
the scaled parameter (see Eqs. (6), (8) and(9)). The dashed lines are the GOP prediction obtained 
from simulations of H(x) = H\cosx + i^sina; with N = 150, using the middle third of the 
spectrum. On the left are calculations in the IBM in its chaotic regime, including 80 states out of 
117 for J = 2 (pluses) and 130 states out of 184 for J = 6 (squares). On the right are calculations 
in the Anderson model in its diffusive regime, including the middle 200 states out of 729. Results 
are given for a cylindrical geometry with W = 4 (circles) and for a toroidal geometry with W = 2 
(crosses). 

FIG. 2. The eigenfunction correlator o(oj,x — x') (see Eq. (10)) as a function of the scaled 
parameter for several values of u measured in units of the mean spacing. The solid lines are the 
analytical results (see Eqs. (13)-(14)), the dashed lines are the GOP simulations, and the circles 
and crosses are the results for the same cases of the Anderson model as in Fig. 1. The 5- functions 
in (10) are regularized with a Lorentzian of width 5 = 0.05 in units of the mean spacing. 
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